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Abstract 

We formulate a thermodynamically consistent equation of state 
(EOS), based on excluded-volume approach, for a hot, dense hadron 
gas (HG). We calculate various thermodynamical quantities of HG 
and various hadron ratios and compare our model results with the 
results of other excluded-volume models and experimental data. We 
also calculate various transport coefficients such as 77/s etc. and 
compare them with other HG model results. Furthermore, we test 
the validity of our model in getting the rapidity spectra of various 
hadrons and the effect of flow on them is investigated by matching 
our predictions with the experimental data. 

1 Introduction 

The search for any unambiguous signal of quark-gluon plasma (QGP) in 
heavy-ion collisions is motivated by properly investigating the behaviour 
and properties of hot and dense HG pQ. Here we formulate a new ther- 
modynamically consistent excluded-volume model where we assign a finite 
hard-core volume to each baryon but mesons in the theory can easily over- 
lap, fuse and interpenetrate into each other. Earlier, we have used this 
model successfully in obtaining the conjectured QCD phase boundary and 
thus determining precisely the location of QCD critical end point [3J [3] . We 
calculate various thermodynamical quantities like number density etc. of 
HG and compare our model results with that of URASiMA event generator 
[4]. We use our freeze-out picture [5] for calculating various hadron ratios 
and compare our results with the experimental data and various excluded- 
volume models. In order to make the discussion complete, we further derive 
r]l s etc. from our model and compare them with other models. Further, we 
extend our model to deduce the rapidity as well as transverse mass spectra 
of hadrons and compare them with the experimental data available in order 
to illustrate the role of flow present in the fluid. Finally, we give summary 
of this work. 



2 Model Descriptions 

Recently, we have proposed a thermodynamically consistent excluded-volume 
model for a hot, dense HG [H |21 Ej. The attractive interaction between 
baryons and mesons is realized by including the baryon and meson reso- 
nances in our model calculation. The repulsive interaction between baryons 
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is modelled via giving an equal and finite size to each baryon. Mesons can 
fuse and interpenetrate into each other so, they are treated as pointlike 
particles. Using quantum statistics, the grand canonical partition function 
for baryons can be written as follows : 



6lT T Jvp Jo ^W+mf [ exp (B^ij + i] 

where gi is the degeneracy factor of ith species of baryons, Ei is the energy 
of the particle (Ei = \/k 2 + mf), V° is the eigen volume assigned to each 
baryon of ith species. Apparently, our approach is more simple in com- 
parison to other thermodynamically consistent excluded- volume approach 
which often involve transcendental expressions and are difficult to solve. 
We determine chemical freeze-out temperature and chemical potential by 
fitting our results with some experimental data and then use them in de- 
termining all other quantities and ratios. 

The rapidity distributions of baryons using thermal source can be writ- 
ten as follows [7J : 
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where is the mass of the ith species. V is the freeze-out volume of 
the system. The rapidity spectra of hadrons with the effect of flow can be 
calculated by using following formula [7] : 
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where rj max is a free parameter related with the longitudinal flow velocity 

Our excluded- volume approach involves hard-core repulsion arising be- 
tween two baryons but mesons do not possess any such repulsion. Wo 
have taken an equal volume V° — 47rr 3 /3 for each baryon with a hard-core 
radius r = 0.8 fm. We have also taken all baryons, mesons and their res- 
onances having masses upto 2 GeV/c 2 in our calculation of HG pressure. 
We have also imposed the condition of strangeness neutrality by considering 
J2i Si(nf — nf ) =0 where Si is the strangeness of ith hadron. 



3 Results and Discussions 

In Fig. 1, we have plotted the variation of number density of hadrons with 
respect to temperature at fixed net baryon density. Our model results show 
a very good agreement with the results of Sasaki 13] except at higher T. 
Figure 2 represents the variation of A/ir~ ratio with ^snn- We find that 
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Figure 1: Variation of total number density of hadrons with respect to 
temperature at constant net baryon density, lines show our model calcu- 
lation and points are the data calculated by Sasaki using URASiMA event 
generator. 




Figure 2: The energy dependence of A/tt~ ratio. Dashed line is the result 
of Cleymans-Suhonen model [8]. Symbols are the experimental data [5]. 
RHIC data are at mid-rapidity. 
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Figure 3: Variation of rj/s with temperature for [1b = in our model and 
a comparison with the results obtained by Gorenstein et al. [9]. 



our model calculation gives much better fit to the experimental data at all 
energies in comparison to Cleymans-Suhonen model [5]. Figure 3 depicts 
the variation of rj/s with respect to temperature as obtained in our model 
for HG [5 having a baryonic hard-core size r = 0.5 fm, and compared the 
results with those of Gorenstein et al. [9]. We find that near the expected 
QCD phase transition temperature (T c = 170 — 180 MeV), -q/s shows a 
lower value in our HG model than the value in other model. In Fig. 4, 
we show the rapidity distributions of 7r + for central Au + Au collisions 
at Jsmn = 200 GeV. Dotted line shows the distribution of tt + due to 
stationary thermal source. Solid line shows the rapidity distributions of 
7r + after the incorporation of longitudinal flow in our thermal model and 
results give a good agreement with the experimental data [7]. 

4 Summary 

In summary, we find that our model results mostly show very close agree- 
ment with those of Sasaki, although the two approaches are completely 
different in nature. Our model calculations for the particle ratios describe 
the experimental data very well. Transport quantities are also successfully 
described by our model. The rapidity distributions which essentially are 
dependent on thermal parameters, also show a systematic behaviour and 
their interpretations most clearly involve the presence of a collective flow 
involved in the final description of the fireball. 
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Figure 4: Rapidity distributions of ir + at v /s]vjv = 200GeV. Dotted line 
shows the rapidity distribution calculated in our thermal model. Solid line 
and dashed line show the results obtained after incorporating longitudinal 
flow in our thermal model. Symbols are the experimental data [7J. 
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